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A treatment of the modal decomposition of the pressure field in a combustor as determined by two pressure time 
history measurements is developed herein. It is applied to a Pratt and Whitney PW4098 engine combustor over a 
range of operating conditions. For modes other than the plane wave the assumption is made that there are distinct 
frequency bands in which the individual modes, including the plane wave mode, overlap such that if circumferential 
mode m and circumferential mode m — 1 are present then circumferential mode m — 2 is not. In the analysis used 
herein at frequencies above the first cutoff mode frequency, only pairs of circumferential modes are individually 
present at each frequency. Consequently, this is a restricted modal analysis. As part of the analysis one specifies mode 
cut-on frequencies. This creates a set of frequencies that each mode spans. One finding was the successful use of the 
same modal span frequencies over a range of operating conditions for this particular engine. This suggests that for 
this case the cut-on frequencies are in proximity at each operating condition. Consequently, the combustion noise 
spectrum related to the circumferential modes might not change much with operating condition. 
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Nomenclature 

amplitude of clockwise pressure wave 

eigenvalue from Eq. ( 14 ) 

amplitude of counterclockwise pressure wave 

resolution bandwidth, Hz, B e = l/T d = 

r/Np = 1 1.71875 Hz 

speed of sound in m/s 

delay time r, s (sometimes expressed as a 

number of samples, D = r48000) 

radial duct mode Eq. (13) 

frequency, Hz 

upper frequency limit, f c = 1/2 At = r/2, Hz 
(24,000 Hz) 

Eq. (17) 

Eq. (18) 

power spectral density 
cross power spectral density 
Bessel function of first and second kind and 
order m 

positive imaginary square root of — 1, \f~l 
cost function 

number of frequencies, f c / A f = Np / 2 (2048) 
duct Mach number 
circumferential mode order 
segment length, number of data points per 
segment (4096) 

number of disjoint (independent) data 

segments/blocks, n d = B e T totai = 234 

number of overlapped segments, with 50% 

overlap n ol = 468 

number of segments 

Fourier transform of p(t) 

pressure signal 
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Superscripts 


conversion constant, 10 (l71 / 2 °)/0.25 

eigenvalue from Eq. ( 15 ) 

contaminating noise at station 0 

contaminating noise at station 9 

sample rate, samples/s (48,000) 

exterior concentric cylindrical duct radius 

interior concentric cylindrical duct radius 

signal at station 1 

signal at station 2 

temperature, K 

record length of segment i 

total record length, s (r^20 s) 

time 

weighting function also known as a window 
function 

axial propagating component of the pressure 
Eq. (16) 

cylindrical coordinates 
analytical coherence threshold 
coherence function 

estimated magnitude squared coherence 
function 

frquency step, 1 /T d , Hz (11.718) 

sampling interval, 1/H1/48000), s 

angular microphone location 

radial duct mode order 

hub-tip ratio, r a / r b 

propagation time delay 

circumferential spinning component of the 

pressure 

phase angle 

frequency in rad/ s ,a> = 2nf 


imaginary part 
running segment index 
real part 
signal x( t) 
signal y(t) 

microphone reference location 


complex conjugate 
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+ = clockwise rotating mode 

— = counterclockwise rotating mode 

= ensemble average 


I. Introduction 

T O UNDERSTAND combustion noise measured in the far field 
of the PW4098 engine it is necessary to investigate and 
understand the pressure field in the annular combustor. The data 
analyzed are from a Pratt and Whitney PW4098 engine. The test was 
conducted as part of phase 2 of the NASA Engine Validation of 
Noise Reduction Concepts (EVNRC) Program. For this test two 
pressure transducers were mounted in the combustor and four far- 
field microphones were used. The coherent combustion noise 
seemed to have a modal pattern. Consequently, a restricted modal 
model was developed and tested to see if it was compatible with the 
available data. This work is reported herein. To be useful, the 
pressure field must be measured and analyzed in a manner consistent 
with the physics of the propagation process. A scheme that 
accomplishes this is discussed in this paper. 

The pressure transducers failed during the test. Consequently, a 
sensor validation analysis was conducted []_]. Results of this study 
were used to select the data evaluated here. The pressure transducer 
signals and the far-field microphone signals are used to obtain the far- 
field coherent combustion noise output power and estimate the 
amount of coherent combustion noise that appears in the far field [2], 
A method for separating correlated noise sources (core noise) and 
uncorrelated noise sources (jet noise) in far-field measurements of 
turbojet engine core noise using multiple microphones has also been 
developed and tested [3], Results from the pressure transducer 
validation study conducted [1] and this study of annular duct modes 
measured in the combustor by pressure transducers are used [2,3] to 
interpret far-field acoustic measurements. 

Combustion noise in a turbofan engine is classified according to 
source mechanism as either direct or indirect. Direct combustion 
noise is attributed to the combustion process. Indirect combustion 
noise is produced when hot spots pass through a turbine or nozzle. 
The existence of hot spots correlated with combustion noise has been 
shown by Miles [4], A review of combustion and core noise is given 
by Mahan [5], More recent developments are discussed by Gliebe 
et al. [6], With either combustion noise source or a combination of 
combustion noise sources, the acoustic pressure field inside an 
annular combustor is governed by the solution of an eigenvalue 
equation. The pressure is assumed to propagate as if it were in an 
infinite duct. Consequently, boundary conditions at the duct ends, 
reflections from the duct ends, and the possibility of standing waves 
are not considered. The pressure propagates in particular bands 
identified by frequency ranges . The pressure pattern in a combustor is 
analyzed by breaking it up into the natural frequency bands of 
pressure propagation in an annular duct. Frequencies in the lowest 
band propagate in a plane wave mode. Higher bands move in 
propagating modes composed of eigenmode solutions to the annular 
duct eigenvalue equation governing wave propagation in an infinite 
annular duct with wall boundary conditions. The modes are 
composed of sinusoids in a circumferential direction and com- 
binations of Bessel functions in a radial direction. 

The physics of acoustic waves in hard circular ducts is discussed 
by Morse [7] (p. 509). The eigenvalue solution for a circular duct 
with flow is discussed by Eversman [8], The physics of pressure 
waves propagating in annular ducts without flow is discussed by 
Tyler and Sofrin [9] in a treatment of axial flow compressor noise. It 
is also presented by Morse [7] as a problem (p. 603). Each higher 
mode propagates only above its own cutoff frequency. It is expected 
that the pressure field in a combustion duct extends over a frequency 
range in which several modes are propagating. The distribution of 
energy in the frequency band for each propagating mode is unknown 
and might provide insight into the physics of the propagation 
process. Consequently, our understanding of the pressure field in an 
annular combustion would be improved if one could determine the 
distribution of energy among duct modes. 


The use of normal mode theory to analyze the propagation of 
noise in ducts to find the optimal microphone positions to estimate 
sound power from sound-pressure measurements was studied by 
Dyer [10], This analysis was soon extended to look at spinning 
acoustic modes generated by fans [9,1 1-15]. This analysis has also 
been used to look at the modal content of noise generated by a jet 
in a pipe [16-19]. An acoustic modal analysis of a YF102 com- 
bustor installed in a ducted test rig was conducted by Karchmer 
[ 20 ] who used the assumption of equal amplitude clockwise and 
counterclockwise spinning modes in his analysis. Equal amplitude 
clockwise and counterclockwise spinning modes result in a 
stationary modal pattern. He found that he only needed the 
circumferential modes 1-6 and the corresponding zeroth order 
radial modes in his reconstruction of the measured data. Acoustic 
modal analysis of the pressure field in the tailpipe of a turbofan 
engine was conducted by Krejsa and Karchmer [21]. and again the 
assumption of equal amplitude clockwise and counterclockwise 
spinning modes was used in the data analysis. 

Before duct mode theory was applied to study core noise other 
techniques were used. Core noise from a Pratt and Whitney JT8D 
was studied by Grande [22] using multiple microphone measure- 
ments to obtain auto- and cross-power spectra in a tailpipe extension 
as well as far-field autopower spectra. Ten internal microphones 
were used to obtain measurements of narrow-band and one-third 
octave band pressure level spectra in an Avco Lycoming YF102 
combustor by Wilson and O’Connell [23-25], Core noise in General 
Electric engines has been studied by Matta et al. [26] to relate 
performance with emissions and noise. Coherence functions and 
transfer functions were used by Moore and Doyle [27] in a core noise 
investigation of the CF6-50 turbofan engine. They used five pressure 
transducers in the combustor, nine other internal sensors, and 15 far- 
field microphones. Moore and Doyle [27] removed the time delays 
between the internal and far-field signals using a cross-correlation 
analysis. Moore and Doyle [27] set the number of averages they used 
to 100 and as a consequence data below coherence function values of 
less than 0.1 are ignored. In this paper, the coherence function is 
plotted on a log amplitude scale because the coherence floor is near 
0.01 . Moore and Doyle [27] tried to identify the source location from 
vectoring of cross-correlation time delays and did not do a modal 
study. However, examination of internal coherence and transfer 
function plots indicate that combustion modes might be present. In 
addition, Moore and Doyle [ 27 ] do state that the appearance of a 
double peak with positive and negative time delays in some of the 
cross-correlation plots suggest that waves are moving circum- 
ferentially around the combustor. 

This study was conducted to understand far-field acoustic 
measurements related to turbofan engine core acoustic noise as part 
of a noise reduction research program. However, related inves- 
tigations of turbofan combustors and annular duct modes are made to 
understand combustion instability. A tutorial on this topic is given by 
Dowling and Stow [28]. An interesting three-dimensional linear 
acoustic analysis of gas turbine combustion instability is given by 
You et al. [29]. The presence of annular duct modes is also a 
fundamental element in the study of fan noise (Envia [30]), fan noise 
radiation from a turbofan engine inlet and aft ducts (Nallasamy [31]), 
and compressor noise (Enghardt et al. [32]). 

The problem of mode propagation in an annular combustor is 
complex due to the presence of high speed flows and temperature 
gradients. The results given here are limited to using circumferential 
modes which resemble plane waves and the plane wave mode to 
replicate the measurements. To analyze the pressure distribution in a 
duct into its modes, it is necessary to be able to measure the 
relationships among the pressures at many different points. 
However, for this test measurements were available at only two 
points. Consequently, in this paper we only investigate the feasibility 
of using a restricted pressure modal model to replicate the single 
measured relationship. 

The investigation discussed will 

1) Show the pressure transducer instrumentation was functioning. 

2) Provide evidence the peaks and dips in the measured 
autospectra and cross-spectra magnitudes are due to duct 
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propagation modes rather than a fluctuating pressure generating 
mechanism. 

3) Identify the energy distribution in the plane wave mode. 

4) Examine the validity of the assumption that one has equal 
amplitude clockwise and counterclockwise spinning modes. 

5) Test the idea that the data can be reconstructed using the first 
few circumferential modes, m =1, m = 2, m = 3, . . ., and the 
corresponding zeroth order radial modes /z = 0. 

6) Determine the energy distribution in each band of frequencies 
composing a propagation mode. 

This paper has two distinct aspects. Because pressure measure- 
ments were only available from two transducers, the modal 
decomposition methodology is restricted to two modes at a time. In 
addition, we must note that acoustic propagation of sound by a 
spinning mode of order m and frequency / will occur only if its 
frequency is greater than a cutoff frequency, / cuto ff . One natural 
outcome of using an appropriate number of microphones is the 
automatic identification of the cutoff frequencies of each mode. 
Consequently, a procedure to estimate the value of the cutoff 
frequency for each mode will be discussed. 


II. Experiment 

The two pressure measurements made in a Pratt and Whitney 
PW4098 combustor will be discussed next. The measurements were 
made in a study of aircraft engine core noise conducted as part of the 
NASA EVNRC Program. The pressure transducer at 127 deg is 
identified as sensor 1 and produces a pressure signal, p 0 (f) while the 
other at 337 deg is identified as sensor 2 and produces a pressure 
signal, pg(t). Sensor angles are measured clockwise from top dead 
center viewed from the rear. The combustion chamber is annular. 
The PW4000 was fitted with 24 injectors [33]. 

The spectral estimate parameters are shown in Table 1 . The signal 
processing algorithms used were written in FORTRAN. They are 
based on subprogram modules developed by Steams and David [34], 
They were modified for this project to provide for a time delay 
selection. In the calculations the segments were overlapped by 50%. 
Results of 10 sets of data measured on the first test day are reported 
herein. This paper shows typical results using four sets of combustor 
measurements. The values of N1 CORR. used in the examples are 
1622 rpm, 1750 rpm, 1900, and 2101 rpm. 

A. Measurement of Autospectra and Cross Spectra 

For discrete time signals the computation procedure is as follows. 
The total record length ^total for received signals pg(t) and p 0 (t) is 
divided into n d disjoint (independent) data segments/blocks, each of 
length T d = T toUi /n d so that there is an ensemble of measurements 
{Pmf 0) and {p 0i (t)} 


which was used as an arbitrary reference 6 = 0 deg location. The 
relative angular location of sensor 2 at 337 deg is then 0 = 210 deg. 

Using notation similar to that used by Piersol [35], the one-sided 
autospectra, G 00 (f), Gggff) and the cross spectrum G 0 g(f) are 
estimated for a given pair of random pressure signals, at a frequency 
/as 

9 "d 

Gooff) = Poifinff) = ... ^/7 i ,(/)/’u,(/) ( 2 ) 

1 d n d 

9 n d 

G m (f) = P 9 ff)P%ff ) = V Pl<J)Pm(f ) (3) 

1 d n d (— 1 


2 nd 

Gooff) = Po(f)PKf) = pt mff) p oiff) ( 4 ) 

1 d n d 

where the”indicates an ensemble average and Pg^f) and P 0l (/) are 
the Fourier transform of the measured pressure signals at the 
reference position and at the angle 6. The * designates the complex 
conjugate. The Fourier transforms used are windowed Fourier 
transforms: 

w(r) = ( 1 - P* ei ff) = f d w{x)p ei {x)e i2 ” fr dx 

V " r 0 ° (5) 

Poi(f)= f ' wfT)p 0i fx)e- J2nfr dx 

Jo 

where w(x) is a weighting function also known as a window 
function. This one is known as a Bartlett (triangular) window. In 
practice these quantities are computed using fast Fourier transforms 
and a nonparametric technique of spectral estimation known as 
periodogram averaging. Several nonparametric techniques of 
spectral estimation are available to be used in evaluating spectra of 
random data (Steams [36], Sec. 15.6 or Hayes [37], Sec. 8.2). These 
techniques differ in the selection of rectangular, Bartlett, Hanning, 
Hamming, or Blackman window functions. In addition, the data 
sequences averaged may be calculated so that they overlap by 
varying amounts. Also, the window functions can be applied to each 
sequence or the whole data history. The Bartlett (triangular) window 
is used in discussions herein and of the aligned and unaligned 
coherence technique by Miles [2] because it provides a simple way to 
convey the idea of a window function. Issues to consider in picking a 
method involve consideration of statistical estimates of bias and 
variance and spectral resolution. The estimate G 0 g(f) is a complex 
number such that 

G ooff) = G^ff) ~ jG lm (f)= I G oe (f) | (6) 


{Pm(t)} = feW + qeft)}i, i=l,...,n d 0 < t < T 

{Poi(t)} = {s l (t) + q 0 (t)} i , i=\,...,n d 0 <t<T 

where qg(t) and q 0 (t) are the contaminating noise at receiver 
location 1. Angular positions 6 are relative to sensor 1 at 227 deg 


where 


I G M ff) I =[GL(/) + GL(/)F 


4>o off) = tan 


~ Glm(/) ~ 

GRe(/)_ 


Table 1 Spectral estimate parameters 


Parameter 

Value 

Segment length, i.e., data points per segment, Np 

4096 

Sample rate, r, samples/s 

48,000 

Segment length, T d = Np/r, s 

0.08533 

Sampling interval, At = 1/r, s 

2.0833 x 10~ 5 

Frequency step, A / = 1 /T d , Hz 

11.718 

Upper frequency limit, f c = l/2Af = r/2, Hz 

24,000 

No. of frequencies, Ly = f c / A / = N p/2 

2048 

Time delay, r = 6323/48. 000, s 

0.1317 

No. of independent samples 

234 

Overlap 

0.50 

Sample length, s 

20 


The signal processing algorithms used were written in FORTRAN 
based on subprogram modules [ 34 ] which were modified for this 
project to be in FORTRAN 90 form. In addition, the processing 
scheme was changed to permit delaying one time sequence with 
respect to another. For these computations, a data segment/block 
record length of T d = 4096/48, 000 = 0.08533 s was used and the 
number of disjoint data segments/blocks was n d = 234. 
Consequently, total record length was T d n d ss 20 s. The analysis 
was done with 50% overlap so n ol za 468. Other spectral estimate 
parameters are shown in Table 1. 

The coherence is given by 


Ylg(f) = 


|G 09 (/)| 2 

Goo(f)Gggff) 


( 8 ) 
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The measured coherence calculated using segments overlapped by 
50% is given by 


Yleif) 


i puf)p a m 2 

£"=i \ p oi(f)\ 2 E"=t \ p m(f)\ 2 


(9) 


Given measurements from two sensors to examine, we can say a 
sensor is providing an incorrect signal if the signals are uncorrelated 
with each other or the correlation changes abruptly in some irregular 
fashion. A procedure developed by Miles [JJ based on the concept of 
aligned and unaligned coherence is used to determine the validity of 
the signals from the pressure sensors. Miles [1] shows that by 
comparing a coherence function calculated using aligned and 
unaligned time histories one may decide if the signals from two 
sensors are uncorrelated or have changed in some irregular fashion. 
The pressure data set examined in this paper was selected using this 
procedure. This procedure compares the normal measured coherence 
Yoe(f)> called the aligned coherence, with one calculated by time 
delaying one signal by an amount D greater than the segment/block 
interval T d . The value of D is chosen so that the two signals are not in 
the same segment/block interval T d , and consequently appear to the 
processing procedure to be uncorrelated. For this study a value of D 
approximately 1.5 times the segment/block interval is used. The 
measured unaligned coherence Foe!/’ D ) is not zero. Even if no tones 
are present, the coherence of two disjoint (incoherent/independent) 
random noise records has a value dependent on the number of 
independent data segments/blocks used to calculate the coherence 
n d . In [2] results were obtained from computer simulation that show 
good agreement with the theoretical estimate of the analytical 
coherence threshold 


yL = i - (i - py (io) 

where we use P = 0.95 and n d = n ot . The coherence threshold y% n is 
discussed by Carter [38,39], Halliday et al. [ 40 1 (p. 247), and 
Brillinger [ 41 ] (p. 317). The coherence threshold Ynn has a value 
which is greater than 95% of the values of the coherence of two 
independent time series calculated using n d disjoint data segments/ 
blocks. For the values shown in Table J_ with a 50% overlap, a total 
record length ^Total of 20 s and a resolution bandwidth B e of 
11.718 Hz, the number of data segments/blocks is n d = 
2A/T Total ss 468. The coherence of random noise for this case is 
Ynn = 0.00639. The concept of aligned and unaligned coherences is 
discussed in more detail in [2], 

In addition, if tones are present in the time signal they will be 
present in all ensembles averaged and appear in the unaligned and 
aligned coherence. In addition to aiding the identification of 
persistent tones, the unaligned coherence provides a reference 
coherence. Small values of the aligned coherence greater than the 
unaligned reference coherence can be easily identified as significant 
because one knows what the uncorrelated coherence looks like. 
However, the major support it provides is in determining if a sensor is 
providing a useful signal. If the aligned and unaligned coherences are 
similar then the two sensors are uncorrelated and one sensor has 
failed. 


B. Autospectra and Cross Spectra 

For the four cases discussed, measured autospectra magnitudes 
and cross-spectra magnitudes and phase angles are shown in Fig. 1 . 
Both a measured and an unwrapped phase angle are shown. 
Although the amplitudes of the autospectra are not remarkable, the 
cross-spectra amplitudes and phases are quite notable. The cross 
spectra have pronounced peaks and dips and the phase differences 
have a tendency to be nearly either 0 deg or multiples of 180 deg. 
These features of the cross spectra are similar to those found by 
Karchmer [20] in his study of YF102 combustor modes. The 
presence of this phase angle feature in the YF102 phase spectra is a 
characteristic of standing circumferential waves and it was found at 
all conditions tested by Karchmer. 


C. Coherence 

The measured aligned and unaligned coherence is shown in Fig. 2 
for the four cases. The coherence functions shown have a lumpy 
appearance which is characteristic of the presence of a modal 
structure and indicates that the acoustic energy is in well-defined 
bands. Although many coherence values are less than 0.1 the 
coherence is still valid. As discussed, the aligned coherence is greater 
than the unaligned coherence indicating the coherence is not due to 
random noise. The random noise coherence value is the unaligned 
coherence value. In addition, the measured unaligned coherence is 
generally below the 95% confidence value for the analytic coherence 
threshold of the random noise which is Ynn = 0.00639 and is shown 
as the black line in Figs. 2a-2d. 

Using a statistical coherence threshold procedure discussed by 
Miles [JJ, each test point was examined for validity and sensor 
degradation. Initially, runs with rpm in the range of 1622-2400 were 
used in the modal analysis. However, eventually runs with 2304 and 
2400 rpm were discarded. This left eight cases. 

D. Phase Angle Standard Deviation 

The coherence function is especially important because the phase 
angle standard deviation can be related to the coherence function. In 
Bendat [42] and in Piersol [35] the random error in the phase 
estimates due to statistical sampling is given in terms of the standard 
deviation of the estimated phase angle 9 O0 by 

where a[6 oe (f)\ is measured in radians, and as used herein n s is 
selected to be the number of segments or blocks used in the spectral 
calculations. For the special case where the term in curly brackets is 
small Eq. (IT ) becomes 

(12) 

\Yoo\V2n s 

where for the unknown coherence Yoe(f)> l* 16 estimated coherence 
Yoeif) f rom Efi- (9) is used. A plot of the standard deviation of the 
phase angle in degrees versus coherence is shown in Fig. 3 for 
n s = n d = 234 and n s = n o( = 468. When the coherence is greater 
than 0.02, Fig. 3 shows that the pressure sensor cross-spectrum phase 
angle standard deviation should be between 20 and 13 deg. This is 
accurate enough if no spinning modes are present so that the 
measured phase angle difference is either 0 or 180 deg. 

III. Analysis 

A. Annular Duct Acoustics 

The annular combustor noise is analyzed in terms of a spinning 
pressure pattern in the gap between two rigid-walled concentric 
cylinders r b < r < r a . However, drawings of combustor hardware 
presented by Hill and Peterson [43], You et al. [29], and Dowling and 
Stow [ 28 ] show this theory is based on an idealized combustor. The 
analysis resembles the fan noise analysis of Tyler and Sofrin [9] 
except that the mean duct flow is included. Propagation of acoustic 
modes in an annular duct with no mean flow is also discussed by 
Converse and Hoffman [44] and Moore [15], The mean axial 
velocity is included in a discussion of axial flow acoustic modes 
propagating in a cylinder by several researchers [ 11 , 16-18 , 45 , 46 ] 
and in an annular duct by Envia [47], 

Tyler and Sofrin [9] resolved the sound field into axial propagating 
modes Z(z), circumferential spinning modes of order m, = 

cos (imcj>) (with m wavelengths in the circumferential dimension), 
and radial duct modes of order /z, JJ) (with /z nodes or zero 

crossings of pressure in the radial direction). 

The radial distribution of pressure in an annular duct is given by 

R(r) = ElZ = J m (a™ "j + (13) 
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B. 






c) 1900 rpm (N1 Corr) d) 2101 rpm (N1 Corr) 

Fig. 1 Autospectrum magnitude for sensor (1) at 127 deg and sensor (2) at 337 deg and cross-spectrum magnitude and phase between sensor (1) and 
sensor (2). 


where J m and Y m are the Bessel functions of the first and second kind, 
respectively, a is the hub-tip ratio, r a /r b , amji is the eigenvalue of the 
radial wave number for the (m, /x, a) mode and it corresponds to k'ff 
used by Tyler and Sofrin [9], Qm/i. is also a function of ( m , /x, a), and 
r is the radius in meters. The mode shape is independent on the duct 
flow Mach number, M dliCt . 

The eigenvalues of the solution are given by the roots of 
J'Jcimlo) Y' m (a ( ^a) _ n 

J m\ u rn ji) 1 m \ u miJ,) 


S£m[i — 


The axial propagating mode is 

iz 


j’M) 

Y’JaZ) 


Z(z) = 




^duct) 


X I M duct £ ± 


k 2 - 


1 -M: 


f>2 \ \ a mii 


duct J " 


(«S) 2 1 1/2 


- 1/2 1 


(15) 


(16) 


where the convention that primes denote differentiation of a single 
argument function with respect to that argument used. Values of 
can be determined from 


where k = u>/ c,a> = 2nf is the frequency in rad/s, / is the frequency 
in Hz, and c is the speed of sound in m/s. From Eq. (16), if k 2 is less 
than ( 1 — M 2 )[(o„^) 2 /rj], the traveling wave solution is attenuated 
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a) 1622 rpm (N1 Corr) 


b) 1750 rpm (N1 Corr) 


0.1 


0.01 


0.001 


0.0001 



-^aligned 

unaligned 

1 


-■^aligned 

unaligned 




y nn =0.00639 



y nn =0.00639 


0.1 


0.01 


0.001 


0.0001 



200 400 600 800 1000 0 200 400 600 800 1000 

FREQUENCY, Hz. FREQUENCY, Hz. 

c) 1900 rpm (N1 Corr) d) 2101 rpm (N1 Corr) 

Fig. 2 Comparison of aligned and unaligned coherence between combustor sensor (1) at 127 deg and combustor sensor (2) at 337 deg. 


with distance. From this relationship the frequency / of a 
propagating disturbance must be greater than / clito ff where 


./critical 


c (a. 


(<0\ 

my.) 


2n 


r b 


(17) 


/cutoff C^duct) — \j ( 1 ~ ^duct (/critical (18) 


Table 2 shows the zero flow eigenvalues (a^), Qmy, and the 
cutoff frequencies, / clItoff = (a\nj L )c /2jtr b for the first five order 
modes and fi = 0 using T = 1619.37 K, a = 0.66, r a = 0.558 m, 
r b = 0.837 m, and c = 806.7 m/s. The topic of cutoff frequencies 
will be discussed further in Sec. HI D where the method used to select 
the modal span frequencies will be presented. 


B. Modal Decomposition Method 

The approach implemented here is based on the mathematical 
modal decomposition technique used by Karchmer [201. The 
analysis starts by assuming that a randomly occurring instantaneous 
pressure disturbance produced by the combustion process sends a 
clockwise and counterclockwise pressure wave moving around the 
annulus. At an angle 6 the two waves add to produce the resulting 
pressure written in terms of annular duct modes: 


M - 1 

Pe(t) = Pj(t) + Pg(t) = e * 

m = 0 

(19) 

where A m and B m are the amplitudes of the clockwise and 
counterclockwise pressure waves associated with the lowest radial 
order of the mth circumferential mode, <f> m is the phase of the with 
mode, j = a/— I, and M is the total number of modes present. Note 
higher order radial modes are neglected. At the reference angle where 

e = o 



Fig. 3 Standard deviation of phase angle as a function of yl„ for 
n s = 234 and 468. 


M - 1 

p 0 (t) = P+(t) + P 0 (t ) = e *" + B m ]e^ (20) 

m = 0 


Table 2 Cutoff frequencies (T = 1619.37 K, r a = 0.558 m, 


r b — 0.837 m,c = 806.7 m/s,M duct = 0) 


Order m 

Node sequence 
no., fi 

(ati) 

(JP) \ c 

J cutoff — 2nr b ’ rtZ 

'xcrtiji 

i 

0 

1.207637 

185.245 

—0.408810 

2 

0 

2.412094 

370.003 

-0.348305 

3 

0 

3.610271 

553.797 

-0.314771 

4 

0 

4.799250 

736.181 

-0.270046 

5 

0 

5.976450 

916.757 

-0.223789 
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Consequently, the cross spectrum PoP* e is 

p o^§ TZzl\A m + BJe»- E^O [A n e-M»-^ + B n e~^^] 

M—l N - 1 

= E H (A - A " + S m SJe^”-^ +nS) 

m=0 n=0 
M—l A'-l 

+ I] Y1 (B - B " + A « 5 '" )e KK -*"- nB) (2D 

m=0 n=0 

After taking the ensemble average the products become correlations. 
The assumption is now made that the modes are independent. 
Consequently, the cross correlations vanish. All terms in the 
averaged version of Eq. (21) for which m ^ n vanish. 

M—l M-_l 

TWe = ]>>?„ + KAn)e ime + + A^,)e- ime ( 22 ) 

m = 0 m=0 

Using the definition = cos (md) + j sin(m#) and noting that 
(A m + B m ) 2 = A 2 m + 2B m A m + B- , (23) 

It follows that Eq. (22) becomes 

PoP§= E»=o ( A m + Bm) 2 cos m6 + j Em=o( A m ~ S m) sin md 

= TWeRe + jP^llm = ( 24 ) 

i^i = md + (m)W 2 (25) 

Equation (24) is the general expression for the complex cross 
spectrum between the pressure measured at the reference position 
and the pressure measured at a given angle 9 relative to the reference 
position. The magnitude is given by Eq. (25). The phase of this cross 
spectrum, then, is 


<Poe(A m -B m ) = tan 1 


(PpPV Im 
(PoPg)Rc 


= tan 1 


Em=o( A l -B 2 ,)sinmd 
T,m=b ( A m + B m ) 2 cos m9 
(26) 


Table 3 Mode span 


Mode, m 

Low frequency, Hz 

High frequency, Hz 

0 

0 

339 

1 

187 

550 

2 

339 

800 

3 

550 

800 


D. Mode Span Frequencies 

The theory of mode propagation and the measurements by 
Kerschen and Johnston indicate that modes cut on at particular 
frequencies [16-19], In addition, examination of the modal 
decomposition of cross-spectrum spectra done by Karchmer [ 20 ] 
shown in Fig. 7 of his paper indicates how the cross-spectrum pattern 
of peaks and dips may be considered to be composed of the 
superposition of cut-on duct modes. The frequency range over which 
each mode dominates also shows up in the cross-spectrum phase 
angle measurements at 180 deg shown in Figs. 4(e) and 5(b) of 
Karchmer’ s paper [20]. Consequently, the second assumption made 
herein is that one can choose cut-on frequencies for each mode used 
to model the system. The cut-on frequencies were selected with 
annular duct acoustics theory in mind. However, the major guidance 
in selecting the cut-on frequencies was the measurements of 
coherence and cross spectra. 

The effective frequency limits for each mode will be called the 
mode span. Examination of the modal decomposition done by 
Karchmer [ 20 ] shown in Fig. 7 of his paper indicates that one may 
easily assume that only the most recent cut-on mode is important. 
Consequently, in this model when mode m is cut on then the 
amplitude of the m — 2 mode is set to zero. At any frequency only the 
amplitude of mode m and mode m — 1 must be found. This will be 
designated using the notation E,„ m -i- Because of this restriction on 
the modes, the term restricted modal model is used herein. 

The selected mode span frequencies are shown in Table 3 and 
Fig. 4. A corresponding chart of cut-on mode frequencies as a 
function of mode number, temperature, and Mach number is shown 
in Fig. 5. The cut-on mode frequencies for the m= 1, m = 2 modes 
are not sensitive to temperature or Mach number perturbations. At a 
particular temperature and Mach number many combinations of 


C. Restricted Modal Model 

Equations ( 24 ) and (26), together with the single cross-spectrum 
measurement from the PW4098 combustor, form the basis of the 
modal analysis feasibility study conducted here. The approach used 
by Karchmer [20] assumed that only the first six modes make a 
significant contribution to the pressure field in the combustor over the 
frequency range of interest in core noise studies. Consequently, to 
determine six unknown values of A m and six unknown values of B m 
one would use measurements of cross-spectrum magnitude and 
phase at 12 positions (0) to create 12 equations so that one could solve 
for 12 unknowns. However, in the study of the YF102 combustion 
modes Karchmer [ 20 ] found that the measured cross-spectrum phase 
was either 0 or 180 deg depending on the sign of the real part of the 
cross spectrum. Consequently, from Eq. (26) the imaginary part of 
the cross spectrum is zero and in evaluating modal constants A m and 
B m he could assume A m = B m . Consequently, he could then infer that 
the sources are such that an instantaneous pressure disturbance 
created in the combustor sends equal amplitude clockwise and 
counterclockwise pressure waves traveling circumferentially around 
the combustor annulus. This creates a standing circumferential 
pressure wave rather than a rotating or spinning circumferential 
pressure wave. For this case, the pressure sensor cross-spectra phase 
angles would have only values of 0 or 180 deg. Using this 
information one only needs measurements of cross-spectrum 
magnitude at six positions ( 8 ) to create six equations so that one can 
solve for six unknown values of A m where A m = B m . Consequently, 
the first assumption made is that for this case the pressure is a 
standing wave and A m = B m . This is done even though some of the 
phase angles measured seem to have intermediate values. 



0 100 200 300 400 500 600 700 800 

FREQUENCY, Hz. 


Fig. 4 Frequency span of modes. 



Fig. 5 Cutoff frequencies as a function of Mach number for T = 1600, 
1650, and 1700 K. 
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temperatures nearby and Mach numbers nearby produce cut-on 
frequencies in the same general neighborhood. Because of the lack of 
knowledge of the exact geometry and operating conditions of the 
combustor, the same mode span table is used for all test conditions. 
Adjustments were made to the table so that the values in the table 
produced reasonable agreement with the theory and the measure- 
ments. The temperature range used was selected based on the 
combustion product mixture leaving the combustor and entering the 
turbine as shown in Fig. 6.22 of Hill and Peterson [43], 

E. Other Assumptions and Solution Procedure 

The final assumption is that the reconstructed model autospectrum 
l/’o/'SI 2 should be less than or equal to the measured autospectrum. 

Using these assumptions we have a measurement of cross- 
spectrum magnitude and phase and two equations and two 
unknowns, A m (/) and A m _, (/), at each frequency. The two 
equations are written in complex notation as 

I G oe (f)\e^eV) = |^P|(A m (/),A m _ 1 (/))|^^W,A m - 1 OT) (27 ) 

where 


PqP*s (A m (f)’A m _ l (/)) = £ 4 Ai cos (mff) (28) 

m,m— 1 


I J2^loo S (me)\ 


M(4n(/M m -i(/))l 2 = V 4 Ai cos (mff) (29) 


7 , „ , , . r f 0 or 360 deg: P 0 Pt > 0 

<Poe(f) — <Poe(A m > A m -i) — | 180 d eg: T\P* < 0 (30) 


I PqPq (A m if ) . A m _i (/)) | 2 = 


' m,m—\ ’ 


(31) 


where a mode m is selected to start at a particular frequency. 

The assumptions made that only modes m and m — 1 are active at a 
frequency, that a correct particular frequency can be selected for a 
mode to cut on, and that the pressure wave is stationary and not 
spinning are not perfectly valid. In addition, the measured cross 
spectrum required as input to the two equations is experimentally 
determined and subject to nominal experimental error and statistical 
uncertainties. Also, the measured autospectrum includes random 
noise in addition to coherent signals from propagating waves. 
Consequently, a solution method that provides an optimum solution 
in a least squares sense without derivatives was used. Algorithms for 
minimization without derivatives are discussed by Brent [48] . The 
search technique used in this study is described by Powell [49] and 
FORTRAN computer code is given by Shapiro [ 50 1 and Kuester 
[51]. The code used was a modified version of the one given by 
Shapiro [50] which was updated to be in a FORTRAN 90 style. The 
cost function used is written in terms of sound pressure level and 
phase angle. We have 

L oe (f) = 101o glo (|G oe (7)| 2 ) (32) 


Cl(f) = [L oe (f)-H oe (m,f)] 2 (37) 


C2(f) = 


[Lxx(f) ■ 


0 : 

7f 00 (m,/)] 2 : 


H w (m,f) < L xx (f ) 
Hoo(m,f) > Lxx(f) 


(38) 


C3(/) 


min [(<j>oe(f) ~ 360) 2 , (flV(/) - 0.0) 2 ] 
min [(4> oe (f) - 360) 2 , ( <M / ) - 0.0) 2 ] 

(ke(f) - 180) 2 


if <p oe (m,f) = 0 

if <p oe (m,f) = 360 
if 4> 0 g(m,f) = 180 


(39) 


C = Cl(f) + 1000 C 2 (f) + C 3 (/) (40) 

where 


l O (171/20) 

Pk ~ 0.25 

and Eq. (38) defines artificial constraint C2 which limits the 
magnitude of the autospectrum calculated with the restricted modal 
model and insures that it is near the measured autospectmm. 

IY. Results 

A. Reconstruction of Measured Data 

One measure of the acceptability of the modal amplitudes is how 
well the coefficients reconstruct the measured cross spectrum using 
Eqs. (28-30). Results are shown in Fig. 6. It is clear that the 
procedure discussed is effective in finding a set of coefficients which 
can construct the measured amplitude and phase spectra. 

A second measure of the acceptability of the modal amplitudes is 
the degree the reconstructed acoustic coherent autospectrum based 
on the coefficients [i.e., using Eq. (31_)] compares with the measured 
autospectmm. Results are shown in Fig. 7. The measured 
autospectmm contains hydrodynamic pressure fluctuations and 
random acoustic noise and should be larger than the reconstructed 
coherent autospectrum. At frequencies less than 180 Hz, where only 
the m = 0 mode exists, the results are quite good as expected. At 
higher frequencies the artificial constraint included in the 
minimization procedure is forcing the reconstmcted autospectrum 
to be near the measured autospectrum. 

B. Modal Content 

To characterize the model coefficients the modes are nomialized 
by the maximum value. The normalized coefficients are shown in 
Fig. 8. The mode normalization coefficients are plotted in Fig. 9. 
Figure 9 also shows linear curve fits to the normalization 
coefficients. 


Lxxif) = min[101og 10 (|G u (/)| 2 ), 101o glo (|G 22 (/)| 2 )] (33) 


= 101og 10 (|R 0 P*(A,„(/),A m _ 1 )| 2 ) + 201o glo (p*) (34) 


H 00 (mJ) = 101og 10 (|R 0 P5(A,„(/),A m _ 1 )| 2 ) + 201o glo (p*) (35) 


4> oe (m,f) = (j>oe{A m (J) , A m _ l (/)) (36) 


C. Normalized Amplitude Coefficients 

The normalized coefficients for all cases were studied. A few are 
shown in Fig. 9 with the mode number identified. The uniformity of 
the plots leads to the idea that one could use the mean value of the 
coefficients to synthesize the cross spectra at a range of angles and 
shaft rotation speeds. The mean was constructed using the first eight 
test conditions. Values for test conditions at 2304 rpm (N1 CORR) 
and 2400 rpm (N 1 Corr) were excluded. Mean relative coefficients 
are shown in Fig. K). These mean values were used to synthesize 
cross spectra at various angles and operating conditions using 
amplitudes calculated from the linear curve fit shown in Fig. 9. 
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a) 1622 rpm (N1 Corr) 


b) 1750 rpm (N1 Corr) 






c) 1900 rpm (N1 Corr) d) 2101 rpm (N1 Corr) 

Fig. 6 Measured and reconstructed cross spectrum. 


However, this paper only show comparisons with the measured data 
shown in Fig. 7. 

V. Discussion 

A treatment of the modal decomposition of the pressure field in a 
PW4098 engine combustor as detennined by two pressure 
measurements was developed. A method using aligned and 
unaligned coherence measurements showed that the pressure 
transducer data presented are valid. The treatment uses the 
procedure of analyzing pressure propagation in the combustor in 
the same manner that pressure propagation in a compressor or fan 
is analyzed. A similar procedure was used by Karchmer [20] to 


analyze YF102 combustor coherence measurements made in a test 
stand and engine. This procedure for the operating conditions 
considered uses an analysis in terms of waves that propagate in 
modes such that for waves other than the plane wave the mode of 
propagation is determined by the condition that the frequency of 
the wave is greater than a cutoff frequency. These cutoff 
frequencies are determined by an acoustic eigenvalue equation, 
wall boundary conditions, and longitudinal boundary conditions 
that correspond to an infinite duct with no reflections and are given 
by Eqs. (17) and (J_8). The analysis used herein provides evidence 
that the peaks and dips in the measured cross-spectra magnitudes 
are due to annular duct acoustic propagation modes. The peaks and 
dips in the measured cross-spectra magnitudes are not related to 
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a. 



b) 1750 rpm (N1 Corr) 



« 

a. 



c) 1900 rpm (N1 Corr) 


d) 2101 rpm (N1 Corr) 

Fig. 7 Measured and reconstructed autospectrum. 
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Fig. 8 Relative amplitudes, A m (f) =B m (f). 


longitudinal instability modes which would be at higher 
frequencies. The source or generating mechanism of the 
circumferential acoustic modes observed could be due to either 
direct or indirect combustion noise. 

Figure 7 shows that the measured autospectra are larger than the 
reconstructed autospectra. This was expected because the 
autospectra include random hydrodynamic noise as well as 
correlated acoustic noise. However, Miles [2] has shown the far- 


field correlated combustion noise is only a small part of the total 
noise measured at low frequencies. In addition, Miles [2] has 
shown a great amount of coherent broadband noise does leave the 
nozzle and is observed in the far field at low frequencies. This 
broadband low-frequency noise may be due to indirect combustion 
noise generated by hot spots as they pass through the turbine as 
discussed by Gliebe et al. [6], pp. 209-223. However, it also might 
indicate that the measured autospectra in the combustor consists of 
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rpm (NI corr) 

Fig. 9 Mode normalization factors and linear curve fits. 


a 



Fig. 10 Mean relative amplitudes for m = 0, m = 1, m = 2, and m = 3. 

random hydrodynamic noise, broadband random acoustic noise, 
and correlated acoustic noise, and the acoustic noise in the far field 
is from the combustor. The concept that fan noise consists of tones 
and broadband acoustic random noise is not new and a treatment of 
broadband noise radiation has been given by Rice [52], The 
method was also applied by Rice to the propagation of indirect 
combustion noise discussed in [6]. Unfortunately, the analysis of 
the two point measurements does not cast much light on the noise 
source. 

The comparison of these results for the PW4098 combustor 
mounted in the engine with similar results for the YF102 combustor 
in an engine and installed in a ducted test rig [20] shows many 
similarities in the pressure field because the pressure field for the 
PW4098 and the YF102 depends on the first few circumferential 
modes, m = 1, m = 2, m = 3, . . ., and the corresponding zeroth 
order radial modes /r = 0. Because the YF102 combustor 
component test rig did not include any turbomachinery this 
fundamental structure of the pressure field seems to be independent 
of propagating pressure fields from either the compressor or turbine 
or difference tones created by interacting waves from the 
compressor and turbine. The fact that the energy distribution is in 
frequency bands with well-defined spans (see Tables 2 and 3) 
related to particular cutoff modes of propagation defined by 
Eqs. ( 17 ) and 08) explains the occurrence of well-defined bands in 
Fig. 2 where the aligned and unaligned coherence functions were 
presented. 

The frequency span of the modes for this case does not resemble 
that in the YF102 because the PW4098 is larger than the YF102 
[20], The energy distribution in the plane wave mode and the first 
few circumferential modes is identified over a range of operating 
conditions. Because of the difference in size the energy 
distribution in the modes for the PW4098 is again different 
because the PW4098 is so much larger than the YF102 [20], The 
model has two drawbacks. The first is the assumption that equal 
amplitude circumferential clockwise and counterclockwise 
pressure waves exist for each circumferential mode. The phase 
angles shown in Fig. 6 would be only 0, 360, or 180 deg if this 
were true. This assumption seems more valid for the YF102 data 
[20], The second drawback is shown in Fig. 7 where the 
reconstructed autospectrum is far below the measured 
autospectrum below 200 Hz where the plane wave propagates. 


Above 200 Hz the difference between the measured and 
reconstructed autospectrum is not very large. 

The successful use of the same modal span frequencies over a 
range of operating conditions for this particular engine suggests that 
the circumferential cut-on modal frequencies are in the same range at 
each operating condition. Consequently, the combustion noise 
spectrum related to the circumferential modes might not change 
much with the operating condition. 

VI. Concluding Remarks 

Turbofan engine acoustic measurements are generally made using 
only far-field microphones for certification to compare noise 
reduction technologies. With this type of instrumentation it is 
difficult to obtain sufficient knowledge about noise sources so that 
advances can be made in noise prediction and noise reduction. This 
paper shows that more information on combustion noise can be made 
available if some additional instrumentation is used. 

The treatment of the modal decomposition of the pressure field in a 
combustor as determined by two pressure measurements developed 
here is in reasonable agreement with the measurements. It was 
applied to data from a Pratt and Whitney PW4098 engine combustor 
over a range of engine operating conditions. The method works well 
for the plane wave mode (below 200 Hz). The method gives 
meaningful results at higher frequencies. The method does show it is 
feasible to use the available measurements to produce a restricted 
modal model. 
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